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1 Introduction 

In one dimensional model of wake field generation by a relativistic electron bunch in a 
cold plasma (the beam-generator is modelled by an infinite layer moving in the direction 
perpendicular to its surface) a number of interesting results are obtained. In particu- 
lar, possibility of acceleration of electrons (positrons) by wake fields is shown with the 
acceleration rate proportional to 7q^ 2 Jl|-@ at n ~ (see also the references cited in 
H]) possibility of generation of the fields exceeding the wave breaking limit by a special 
combination of bunches (generator- invertor-dumper) is shown in ||, etc. 

However the obtained data require some refinement on the subject of finite transverse 
sizes of bunches, since the bunches of existing and designing accelerators and storage 
rings mainly have longitudinal geometry, or some times comparable longitudinal and 
transverse sizes. 

In present paper master system of equations for two-dimensional finite charged rela- 
tivistic bunch (of electrons or positrons) moving in cold, collisionless, vortex free, station- 
ary plasma with immobile ions are formulated (see also 0-||). The system of nonlinear 
equations in partial derivatives is then reduced to the form convinient for further numer- 
ical analysis. 



To check the validity of obtained system of equations two limiting cases are con- 
sidered -for infinite wide (transversal) bunch and for infinite long (longitudinal) bunch. 
Both cases admit exact analytical solutions, obtained formely in |Q and [|J. Consid- 
ered limiting cases of the general system as one should expect, coincide with the main 
equations obtained in f2|, [|J. The analysis of solutions of mentioned limiting cases from 
viewpoint of development of numerical algorithm for solution of the main general system 
of equations is also presented. Ranges of parameters values are determined, where non- 
stable or nonphysical (n e < 0) solutions take place, when doing numerical calculations, 
in particular, at joining the solutions in space regions occupied by a limited bunch and 
that free from it, removal of abovementioned limitations on plasma properties, probably 
will promote to the clearing up the nature of these nonstable or nonphysical states. 

The third obvious check of the validity of the numerical calculations is the passage 
to the linear limit, which is considered in a number of papers (see e.g. [|J, 0, || and 
references therein). 

The work is presented in the three parts. Present paper (part I) is devoted to formu- 
lation of the system of the master equations of the problem and analytical considerations 
of the limiting cases, which have a exact solutions [flJ-JJ (transversal and longitudinal 
flat bunches, linear approximation). Presented consideration gives some hints to the 
construction of proper algorithm for numerical calculations. 

The second part (II) will be devoted to the construction of the algorithm for numerical 
calculations of the system of basic equations of the problem for the general case of the 
bunch of finite transversal and longitudinal dimensions. 

The third part (III) will present the results of numerical calculations, their analysis 
and comparison with the existing analytical results. 
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2 The Master Equations 



The equation of motion of relativistic plasma electrons in electric E and magnetic H 
fields is written as: 

f +(v^)p=-eE-l[vxH], (1) 

where p is momentum, v is the velocity of plasma electrons in lab system, the electron 
charge is set — e. Introducing scalar (p and vector A potentials for electromagnetic fields 
through formulas 

1 dA 

E = -gradip - -—, H = rot A (2) 
one can represent the equation (|I|) in the form: 

d 6 — * 6 — * 

— (p A) = egradip — mc 2 grad'y — [rot{p A) x v\, (3) 

where 7 is the Lorentz-factor of plasma electrons. 

One can note, that the eq. @ has a partial solution 

6 -» 

p=-A + gradx, (4) 

mc 7 = eip- (5) 

where % ls an arbitrary calibrating function. Since from @) follows rot(p — ~ c A) = 0, 
this solution corresponds to the vortex free motion of the plasma electrons considered 
in M, M, [Q. The eqs. (fflBl) along with Maxwell's equations for the field potentials, 
with charges and currents corresponding to the motion in plasma bunch charges (which 
is supposed to be given- rigid bunch approximation) form a complete system of equations 
of the cold vortex free hydrodynamic plasma model, with immobile ions. 
Later on we will work with dimensionless potentials 

a = A/(mc 2 /e), (6) 

/ = ^/(mc 2 /e) (7) 
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and we will normalize the density of plasma electrons and bunch charges on plasma ions 
density no- Introduce also plasma frequency uj p and wavelength A p 

uj p = (47re 2 n /m) 1//2 , X p = c/uj p (8) 

we will use dimensionless coordinates (in units of X p ) and dimensionless time (in units of 
Up 1 ); 4-momenta of plasma electrons represent in the form (e,pc) = mc 2 (7 , ^7) , Lorentz 
factor and velocity of bunch charges denote as 70, $oC. 

Let's choose a calibrating function x = 0. As a result we obtain a full system of 
equations for the considered plasma-electron bunch model: 

7 = /, (9) 

h = 3, (10) 
d fdf \ 

Uf+ dt I ~5t +dw3 j = ~ 1 + n e + rib, (11) 

Da — grad + diva^j = n e /3 + nb(3 , (12) 

where n e is the plasma electrons density, n& is the bunch charges density (positive for 
electrons and negative for bunches consisting of positively charged particles). 

As one can see from relations (|9]) and (|T0| ) components of 4-potential are connected 
on mass surface of plasma electrons: 

f 2 = l + a 2 y + a 2 z (13) 



This follows from the selection of calibration function x — 0; i.e. (13) is a calibrating 
condition on potentials. Calibration of potentials (JT3h we will call "energetic". 
From eqs. (|TT| ) and (12) also follows continuity equation for plasma electrons: 



div(n e (3) + -j£ = (14) 

We are interested in steady state wake field solutions of eqs. (|9l-|l"2"D, i.e. only on those in 
which longitudinal dependence of all variables is determined by the longitudinal position 
of point of observation relative to the bunch. Assuming that the bunch propagates along 



the axis z, suppose all the physical variables depending on combination z — fat, which we 
also will denote as z. Let us consider the case when all the physical variables depend on 
transverse coordinate y only (flat bunch with horizontal dimensions much more greater 
than vertical one). Then we obtain three scalar eqs. for three values /, a y , a z and unknown 
plasma electrons density function n e , obeying the continuity eq. fll4]) : 



d 2 f d 2 f d 2 a 

^ + a?-*a? = - 1+ "' + "" (15) 

aS^-*>£ (16) 

d 2 f d 2 a y d 2 a z 2 d 2 a z a z 

(3 °d^-d^y- + W (3 °^ = ne 7 + ^ {7) 



It is convenient to introduce the new unknown functions 



v = f- (3 a z , (18) 
u = Pof - a z 



Then eqs (p~5| - [l~TD convert to the form: 



U^ 1 -®*)--**'* 1 -®^ (19) 

% + (1 - ®>% = - 1 + n ' (1 - ;3 °> + nA1 - (20) 

Denote the combination involved in both eqs. fll9|) and (]21|) by 

>=%Hl-®%, (22) 



one can symmetrize derivatives in ( |19| - |2T| ) introducing "tilded" argument and variables: 

z = 270, v = vj , u = uio, P> = mo (23) 

Instead of n e introduce a modified electron density iV by the formula: 

N = (24) 
7o(w - PqU) 
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The system of equations for variables a y , u, v, jl, N describing considered plasma elec- 
tron bunch model, is written as 

+ || = 7o(-l + Vl-/3 o 2 )) + ^ (25) 

du da v _ , . 

3-y + li = » (26) 

^ = -Polo + Nu (27) 

f = Na y (28) 

1 + a 2 y + u 2 = v 2 (29) 

The last equation is the condition of "energetic" calibration of potentials ([13]). 
One can easily see that in this representation the continuity equation (|14D is resulted 
from (|27|), (|28| ) and has the form 

^{Na y ) - ^(Nu) = (30) 



At numeric modelling of the system of equations (|25|-p9|) one of eqs (0), (|28|) can be 
replaced by the eq. (|3"0|). 

One can reduce eqs ( p5H29]) to eqs. for a y , v, jl. To do so, one can express N e.g. from 
( p7| ) and define u from (|29|). Eventually we obtain the system of equations 
da,, r-z da.,, „dv 



a y^-^ 2 - 1 - a l-^i = %-^ 2 - 1 ' a l-~^ ( 31 ) 

du r— da , . 

<9 2 v <9 2 -u , v f n dil\ /00 , 

V + 5F = w ° " "<° + 7F=T^J + d-y) (33) 

Below we will use bunch models symmetrical with respect to y axis. In this case one 
can easily note, that functions a y and fx are odd functions on y, v is an even one. 

Also note, that nondisturbed plasma (no charged particles bunches and wake fields) 
corresponds to the following values of variables: 

v = 7o, u = Polo, N = 1, a y = jt = (34) 



It is also interesting to write the equations in the form when the calibrating condition 
(p9|)is an integral of motion of the differential equations. This can be done by involving the 
plasma electrons density among the unknown functions. Monitoring of this parameter is 
useful also from the physical viewpoint. Let's write these eqs. for deviations of functions 
u, v and N from the vacuum values: 



u = Pi + X, 



(35) 



v — 7 + A, 
N = 1 + v. 



(36) 
(37) 



Using eqs. (p5j - p8[) , as well as the continuity eq. (|29"D for the values a y , X , A, v one can 
obtain the following eqs.: 



i/)(Ao 



y a y) 



2 dv da y dv dx dv dx 
dy dy dz dy dy dz ' 



(38) 



(1 + v)(A X -x) = K 1 + v)(Polo + X) ~ (Polo + + a VQ~^~ ~ 

^ dv dx _|_ dv da y ^ dv da y 

dz dz dz dy dy dz 

n b 



AA - A = ^ + 1/(70 + A) 
7o 



,(Pv 
' dy 1 



d 2 v 

- 2a y (/5 7o + xYt^qZ + (Polo + x) 



,(Pv 

dz 2 
, dv da,„ 



-(1 + u)( l0 + A)^ + (1 + v)(-v + 7o (A - p oX )) + (Polo + x)(^^ + 
lo dy dz 

| duday _ 2 dv_dx^ + a ,dv_d]C + dv dx _ ^dv da y ^ 

dz dy dz dz y dy dz dz dy dy dy 



-(1 + v) 

The calibrating condition is written as: 



dy) \dz) \dy) \dz J \ dy J \ dz 



(39) 



(40) 
(41) 



27o(A-/?oX) + A 2 -<- X 2 = 
7 



(42) 



In spite of rather unwieldy form of eqs. ( |38| - |4T| ) they have a number of advantages for 
numerical calculations. Firstly, the relationship ( f42l) follows from eqs. (|38|-[|I]) and it can 
used for verification of calculations. Secondly, from eqs. (|38l - |4T|) one can easily come to 
the linear case corresponding to the bunches with low density. In so doing the linearized 
eqs. are written as: 

Aa y - a y = /3o7o7ttt^ (43) 



Ax - X = A)7o(^ ~ ~^r) (44) 

AA - A = — + 7o z/ (45) 
7o 

$707^ = ~n b -y (46) 



The calibration relationship 



A - PoX = (47) 

can be easily obtained from eqs. ([H]4J6|). 

The linearized eqs. are also a convinient testing area for numerical methods in going 
to the nonlinear eqs. 



3 Limiting Cases of the Master Equations 



Let us consider two limiting cases of equations (|25|-f27j): previously considered case of 
bunch with infinite transverse sizes and case of bunch with finite transverse and infinite 
longitudinal sizes. One can note, that in both cases from the formulation of the problem 
follows, that there is no plasma electrons flow along the axis y, i.e (3 y = and 



a y = 0. 



(48) 



Taking into account (H) one can reduce (|3l| - |33D to 



dv r~ 
v- — yv 2 -lu — 0, 
dy 



oz 



(49) 
(50) 



a? + 9F = "«> - 70 + TJ=7 («™ + % j ' (51) 



In the case of infinite transverse size there are no dependence of all variables on y. 
This case is described by a single equation for v: 

|| = W - 7o + Po 7q -^= (52) 

Analoguous equation was obtained previously in |fj] and was a subject on analysis in 
a series of papers e.g. ||. Note, that condition (|) of absence of vortex 

rot(p - -A) = 

is fulfilled automatically. 

If substitute v and z by v and z from (pl|), then eq. ( [52] ) coincides with eq. @ in [||], 
having an integral 

e = V 2 + 7 > - - 7" 2 ) V2 ] - n 6 i;, (53) 

which allows to consider equation of motion of plasma electrons and Maxwell's eqs 
(Coulomb's Law) as an equation of motion of a point with a unit mass, coordinate 
v, velocity v', moving in potential 

U = 7 > - Po(v 2 - 7" 2 ) V2 ] - n b v. (54) 
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On Fig. 1 U is represented as function of v for different values n b and 70 = 10. 
Negative values of n b correspond to bunches of positively charged particles. Boundary 
values v = 1, U' = 0, e = 1 — n b are reached at front of bunch z = d. One can see from Fig. 
1, that for rib < the motion is always periodical, it is also periodical for < n b < 1/2 
and is nonper iodic for n b > 1/2. 
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Fig. 2c. Phase-plane portrait for n b =10, y =10. Separatrix: 5=9. 



On Fig. 2 phase space portraits of system cold plasma -one dimensional flat bunch 
of length d with infinite transverse sizes, rib = — 10; 0.4; 10,7 = 10; boundary values are 
not fixed (different e) are presented. One can use these trajectories, in particular, to 
form solutions for combined bunches moving in plasma and generating no wake fields M. 
Solution for a bunch with given d, rib, A) must be joined with wave solution for a free 
plasma beyond the bunch (wake waves with rib = 0, phase velocity of waves coincides 
with velocity of the bunch v p h = v o [|TJ, @). In doing so no additional limitations arise if 
rib < 0, or < rib < T+8~' ^ or non P er i°dical solutions inside bunch, when rib > y+" 



may become negative; it remains positive for values of dimensionless pulses p(0) = p ez /mc 
of plasma electrons on the rear edge of the bunch —Polo < p(0) < 0, i.e. for certain 
bunch lengths < d < d cr ,p(0) is limited from above. Negative n e usually interpreted 
a consequence of violation of steady state condition(see e.g. @, 0); in particular, it is 
possible that the wake field will break behind the bunch if the wave amplitude is large 
enough imidiately behind the bunch. 

In the other (longitudinal) limiting case there is no dependence of all variables on 
z. Physically such conditions can be realized for long enough bunches far from the 
bunche's head. The linearized problem for semi-infinite bunch in cylindrical geometry 
was considered in fl2| . A tensor of dielectric permeability was introduced to connect the 
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plasme electrons current density with electric field. The essential contribution was defined 
by concomitant fields. The transition fields arising at injection decrease exponentially. 
In the case when the bunch transverse size is much smaller than the plasma wave length, 
a return current, concentrated in the bunch transverse section, arizes. In the other 
limiting case, the transverse section of the return current is much greater than the bunch 
transverse section. 

From (f49|-|5"l|) we obtain 



d 2 v _ x v d 2 



, ?Vfa -70+ (PolO + TT^yv 2 -l) (55) 

oy J v 2 _ i ay 1 



It is interesting to note, that eq. (|5"5|) (the longitudinal limiting case) differs from (|52j) 
(the transverse limiting case) by the term d 2 \/v 2 — 1/dy 2 in brackets in the right side of 
eq. (H). 

By substitution 

Vv 2 - 1 = shR, v = chR (56) 
we obtain the following equation for R 

R" = - lQ )shR - (3 0l0 chR (57) 

Prime in ( f£2"D means derivative with respect to y. Note that by substitution 

R = ln l0 {l + [3 )- X (58) 
one can obtain following equation for \: 

X" = n b p chx + (1 - n h )shx (59) 

This equation is obtained in work j3J, dedicated to the transportation of longitudinally 
infinite bunches of charged particles with current, exceeding the Alfven limit through a 
cold vortexless plasma. 

The first integral of ([5*9]) when = const is 



E = -{R!f + (n 6 7 -1 - -? )chR + M shR = const (60) 
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Thus in this case also the problem is reduced to the study of massive (unit mass) 
point motion in potential field U(R): 



U(R) = -flposhR - {n b % 2 - l)chR\ (61) 

Form of function U(R) essentially depends on the regions of parameter values rib and (3$. 
When 

n b > j^-r, (62) 

the function U(R) has a single minimum and form the potential well. 
If 

1 . . 

i + Po 

the function U has one maximum and looks like a "hill". 
In region of values 

< n b< T—Ti ( 64 ) 



i + A) i-/V 

the function has no extrema and has a characteristic bend, and varies from — oo to +oo. 




In the fig. 3 region of parameters (rib, Po) and some potential curves for (3q = 0.85 
and a few values of rib are represented. 
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Phase trajectories in the plane (R,R f ), corresponding to parameters (rib, Po) from 
regions I (n b > 1/(1 - #,)), II (1/(1 - #,) > n b > 1/(1 + /3b)), HI (1/(1 + #,) > ra 6 ) are 
presented in figs. 4-6. 





Fig. 4. Phase trajectories lor parameters tram region I are shown (n b =7, pV0.85) 
Curves tor energy values from E=12.78 down to 0.78 are represented; E=0.77 
corresponds to the immobile point of phase trajectories. 



Fig. 5. Phase curves for region II (n b =0.7, hVO.85). Phase curves corresponding to 
the values of energy E=6, 4, 2, , -2, -4, -6 are shown. There are no immobile 
points lor these phase curves. 




Fig. 6. Phase-plane portrait for region III, n b =0.4 (B =0.85). The immobile point is 
defined by the value R=1 .898, the energy E=-0.494. Other phase curves 
correspond to the shift o( energy by values ±2, ±4, ±6. 



For the region I the phase trajectories are closed and there is a single stable immobile 
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point: 



R' = 0, Rq = arcth ^° _ 2 (65) 

In the region II all phase trajectories are not closed, there are no singular points, in 



the region III only hyperbolic motion around the point (65) can be realized; the point 
(|65|) in this case becomes unstable one. 

Location of Rq for different values of parameters rib are shown in fig. 7. 
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Fig. 7. Location of immobile point of phase trajectories on the R axis depending 
on n b . Region of allowable values thR is defined by the strip (-1 , +1 ). 



Since finally we will be interested in restricted in transverse direction bunches (rib = 
at \y\ > b,2b is the width of flat bunch), the case rib=0 requires an additional consid- 
eration. Passing from the region occupied by bunch (\y\ < b) to the region without it 
needs the physical solutions, corresponding to the null value of all physical quantities 
(electromagnetic fields, pulses of plasma electrons) at y — >• ±oo. These are the solutions 
which must be joined with the solutions inside the bunch, where rib 7^ 0. For hyperbolic 
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phase trajectories (corresponding to the case rib = 0) the separatrices passing through 
the unstable (hyperbolic) immobile point are the only suitable trajectories. Passing along 
two branches of separatrices requires infinite long "time-development" on parameter y. 
Eqs. for the branches of separatrices are: 



R' = ± 



1 ~ Po \ V4 e R/2 _ ( 1 + Po \ V4 e -R/2 



(66) 



1 + PoJ V-Po / 
Construction of solutions, corresponding to the bunches finite with respect to y, can 
be done in the following way. Branches of separatrix rib — 0, realizing the motion in y from 
the side surfaces of the bunch to the infinity are to be superimposed on phase trajectories 
corresponding to the position of ((3o, rib) in one of the regions I, II, III. Separatrix branches 
intersect many phase trajectories-different choices of intersection regions correspond to 
the different widths of bunches with the same values of parameters (rib,Po)- Figs. 8 and 
9 represent the process for construction of solutions in regions II and III. 




R 



Fig. 8. Formation ol physical solutions for limited along y bunches; MM, p =0.85 Fig. 9. Formation of physical solutions for limited bunches (n b =0.7, (i =0.85, II 

Branches ol separatrix for n b =0 are shown, which realize transition to ±» along y axis. region) analogous to Fig. 7. Arrows show the limiting combined phase trajectory 

One of combined trajectories is shown, the interval of phase trajectory between corresponding to the plasma electrons density n e =0 on the R axis, 
branches of separatrix correspond to the solution inside the bunch. 



Note an important circumstance. There are no transverse motion of plasma electrons 
in adopted model and adjacent layers of plasma slide one with respect to other (because 
we neglect the plasma viscosity). In such a formulation the continuity equation is satisfied 
automatically and from the mathematical viewpoint there are solutions with any values 
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of n e , including negative ones, found, in particular in 0, ||]. 

However, clearly such solutions are not physical and are to be discarded. In [|J an 
illegitimate, (by our opinion), attempt is done to interpret the regions with negative n e 
as channels from which all plasma electrons are displaced. 



Using the formula fl24|) and relation (pOf ) one can obtain the following expression for 
plasma electrons density: 

n e = 7o c/lR(1 - (3 thR)(-E + ^(R'f). (67) 

From (|67|) follows that n e can alter its sign, and it can take place only if E > 0. 
Indeed, as it is seen from construction of physical solutions by means of joining of phase 
trajectories for rib 7^ with branches of separatrix at n& = 0, such a combined phase 
trajectory always intersects the R axis (R' = 0). Thus, for positive E inside the bunch 
n e < 0. It is clear, that there are no physical solutions when rib > 1/(1 — Po) (region I 
of values (n&, A))), since the minimum values of corresponding potential pits are positive. 
In contrary, in the region III (rib < 1/(1 + /3q)) the condition E < does not impose any 
additional limitations, since the regions for joining of solutions correspond to E < E Q , 
where Eq < is the maximum value of potential U(R). 

Lastly, in the II region (1/(1 — (3q) > rib > 1/(1 + A))) intervals of phase trajectories 
necessary for joining with branches of the immobile point of phase trajectories correspond 
to the energy values E > 0, when rib > 1- Therefore at 1/(1 — /3 ) > n b > 1 inside bunch 
always n e < 0. Such solutions as it was mentioned above are not physical by our opinion. 
In the region 1 > n& > 1/(1 + @o) the condition n e > is not satisfied for a limited set of 
values of energy only (rib — 1 < E < 0). The maximum allowed bunch width corresponds 
to the "time" which is necessary for passing the section of phase trajectory at E = 
between the branches of separatrix, where rib — 0. 

In fig. 9 arrows indicate the limiting combined phase trajectory corresponding to the 
plasma electrons density n e = in the middle of bunch. Phase trajectories right to the 
limiting one correspond to the nonphysical (n e < 0) solutions. 

Complete sweeping of plasma electrons by the bunch, i.e. formation of a channel in 
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plasma take place at (-R') 2 = \E — 0- I n particular, this condition, as it was mentioned, 
is satisfied in the middle of bunch (R' = 0) at E = 0, i.e. at R = arcth((3 /(l — rib^Q 2 )) 
in an immobile point of the region I, (which is unstable for region III). 

Thus, summarising the contents of the presented consideration one can note, that for 
a transversely infinite bunch at rib < 1/(1 + (3q) physical solutions exist for bunches of 
any width. In the case 1/(1 + (3q) < rib < 1 physical solutions exist for bunches of limited 
width only. When rib > 1 in the middle of bunch n e < and we discard such a solutions 
as nonphysical ones for considered formulation of the problem (collisionless, stationary, 
vortexless, cold plasma). 

Very likely that consideration of the problems which takes into account violation 
of stationary condition, thermal motion of plasma electrons and ions, plasma viscosity 
possibly will allow to clear out the physical nature of these nonstable or "nonphysical" 
solutions. 

It is necessary to remember that states with n e < are obtained in very long or very 
wide bunches. The finite dimensions of real bunches can change the domain of existance 
of these states or even eliminate them completely. Linear approximation to the limiting 
cases of formulated problem, discussed particularly in M, M, M, is valid, when bunch 
density is n& <C 1, and plasma electrons are nonrelativistic \j3\ <C 1. The last assumption 
provide that the condition rot(p— -A) = for the absence of the vortexes in plasma flow 
is fulfilled automaticaly, due to nonrelativistic equation of motion of plasma electrons 
and Faraday's law ||. Condition n& <C 1 allows to seek the plasma electron density as 
n e = 1 + n' e , \n' e \ <C 1 and linearize the continuity eq. (|14|) and then, using nonrelativistic 
equation of plasma electron motion and Coulomb law, find a solution for n' e . According 
to the condition (3 <C 1, / in @ in linear approximation is / = 1 + /', |/'| <C 1, a ~ (3 and 
eqs. ([LTD , dl2|) with the given right side have exact analytical solution, which can served 
as a test function for subsequent limiting case in computer simulation of the general 
problem. It is necessary to take into account that condition /3 <C 1, can be fulfiled only 
for short enough bunches. 

Obtained results of presented analytical consideration of general system of equations 
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for two-dimensional limited bunch of charged particles moving in a cold, collisionless, 
vortexless, stationary cold plasma, as well as the limiting cases corresponding to the 
infinitely wide (transverse) and infinitely long (longitudinal) bunches and linear approx- 
imation, valid for nj, ■€ < 1, must serve as an analytical bases for formulation of 
an algorithm for numerical solution of the problem-obtaining of wake fields and focusing 
forces for bunches of an arbitrary longitudinal and transverse sizes. 
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